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Abstract 

By employing conformal mappings, it is possible to express the solution of certain boundary 
value problems for the Laplace equation in terms of a single integral involving the given boundary 
data. We show that such explicit formulae can be used to obtain novel identities for special func- 
tions. A convenient tool for deriving this type of identities is the so-called global relation, which 
has appeared recently in a wide range of boundary value problems. As a concrete application, 
we analyze the Neumann boundary value problem for the Laplace equation in the exterior of the 
so-called Hankel contour, which is the contour that appears in the definition of both the gamma 
and the Riemann zeta functions. By utilizing the explicit solution of this problem, we derive a 
plethora of novel identities involving the hypergeometric function. 



1 Introduction 

Simple boundary value problems for the Laplace equation in two dimensions can be solved by the 
powerful technique of conformal mappings. However, the usual implementation of the conformal 
mappings presented in most text books is limited, because it applies only to the case that the given 
boundary data is piecewise constant. In section[2]we present a technique which can be applied to the 
general case of arbitrary boundary data (see also [24]). By employing this technique, we construct 
in section [2] a solution of the Neumann boundary value of the Laplace equation in the domain D, 
which is the exterior of the so-called Hankel contour H (see figure [Ll]) defined by 

H = {z = re~ ia , a < r < oo} U {z = ae~ %e , -a<6*<a}U{z = re ia , a < r < oo}, (1.1) 

where 

< a < 2ir, \a\ < ir. 



A solution of the above Neumann boundary value problem is expressed in Theorem 2.2 in terms of a 
single integral involving the three functions {g+(r), g(<p), g_(r)} defining the Neumann data. Using 
this explicit solution, we derive in section [3] four integral identities. A particular case of these four 
identities can be derived by employing the explicit solutions 

q{r,9) =Re(z k ), q(r, 9) = lm(z k ), z = re i6 , r > 0, 9 G R, Rek>0. (1.2) 



In the general case, these identities are derived using the so-called global relation. We recall that 
this relation plays a crucial role for the contruction of the solution of a wide class of initial-boundary 




Figure 1.1: The Hankel contour. 



value problems; for evolution and elliptic PDEs respectively, see for example [61 \TT \ \12 \ \17 \ [T9ll2"5[ I2"B] 
and [21 O El [9l HH [211 HEH EI]. Here, instead of employing the global relation to solve the given 
boundary value problem, having already constructed the solution via conformal mappings, we use 
the global relation to obtain the integral identities mentioned above. 

Regarding the global relation, we note that in addition to its basic role for the analytical solution 
of a large class of PDEs [13l[[H[T5], it nas a ^ so been utilized in the following contexts: (a) It yields a 
novel non-local formulation of the classical problem of water waves [H HI fTU], [T%1 123|: (b) it provides 
a useful approach to Hele-Shaw type problems |7j; (c) it gives rise to novel numerical techniques for 
elliptic PDEs in the interior of a convex polygon [2U1 l22l [271 1281 [29] . 

2 A Neumann Boundary Value Problem for the Laplace Equation 

The Dirichlet and Neumann boundary value problems for the Laplace equation in the exterior of 
the Hankel contour can be solved via conformal mappings. In this respect we note that the usual 
implementation of the conformal mappings presented in most text books fails here, because it only 
applies to the case of piecewise constant Dirichlet or Neumann data. For the case of arbitrary 
Neumann data the following result is useful. 

Proposition 2.1. Let q(r,6) satisfy a Neumann boundary value problem, for the Laplace equation. 
Let uj(z) be the conformal mapping of the relevant domain to the upper half of the complex u-plane, 
where 



oj(z) = x + iy 



z = re 



x,y 



e 



r > 0. 



(2.1) 



Then, 




(2.2) 




(2.3) 



q y (x, 0)dx 



1 



qgdr — rq r d6 



(2.4) 



r 



2 



Proof. Observe that 

1y dx \ y =o = - qu,du)) y=0 . 

Indeed, the RHS of this equation equals 

i 



(2.5) 



[(Qx ~ iq y )dx - (q x + iq y )dx] 



y=o> 



which equals the LHS of (2.5). 

The conformal mapping implies 

Furthermore, using z = re l9 , we find 



,du) 



rdZ. 



ydZ 



-i6 



2 

[q r dr + 



i 

—t 
r 



+ 



e ie (dr + ird9) 



rq r dO qgdr 

r 



Hence, 



q w du - q^dco = i 



rq r d9 qgdr 

r 



and equation (2.5) becomes (2.4). 

Equation (2.2) is the well known Poisson formula for the Neumann problem. 

Theorem 2.2. Let the domain D be defined by 

IT 

D = {a < r < oo, —a < 6 < a, }, < a < 2ir, — < a < ir. 



□ 



(2.6) 



Let q(r, 9) solve the following Neumann boundary value problem for the Laplace equation in the 
domain D: 

qe(r,-a) = g-(r), a < r < oo, (2.7a) 

q r (a, 9) = g(9), -a < 9 < a, (2.7b) 

qo(r,a) = g+(r), a < r < oo, (2.7c) 

where the functions g±(r) and g{9) have appropriate smoothness and decay. A solution of this 
boundary value problem is given by 

1 2 



q(r,9) 



1 

2^ 



+ 



+ 



1 

2^ 



2tt 



30 9+(p) 
P 

9-{p) 



In 



r\2a 

a ./ 



r\~2a 



cos 



tt9\ 



In 



P 



a 

r\ 2i 
a / 



a 

r\ ~i 



2 2 V2a/ 



1 2 



COS 



2a/ 



+ 



/7T(/?\ i?(r) 



sin — 
V2a/ 



sm 



2a7 



where R(p) is defined by 



R{p) 



a/ Va/ 



dp 
dp 

(2.8a) 
(2.8b) 



and principal value integrals are assumed if needed. 
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Proof. It is straightforward to verify that the function ui(z) defined by 



(2.9) 



maps 



the domain D defined in (2.6) to the upper half complex w-plane. The points {ae ia ,ae ta ,a} 



are mapped to the points {—1,1,0} respectively. Also 



x 



1 . (*0\ 

- - sin — 

2 \2aJ 



r\... / r 
a J \a 



1 (tt6\ 
y= 2 C ° S UJ 



a J 



r \ - — 
a 



(2.10) 



The condition y = implies either = ±a or r = a. Hence, 

— aq r (a, if) dip 

gy(£,0)de= 1 

-qe{p,±a)dp 



Thus, equation (2.2) becomes 

1 



dp 



a dp 1 f°° 
( r ,d) =77- / lnF(/3,a)g e (p,a) h— / hiF(p, -a)qg(p, -a 



1 f a 

— J hxF(a,<p)[- aq r (a,cp)]d(p, 



where 



iW)= [£(p,^)-x(r,0)] 2 + y 2 (r,0) 

and for convenience of notation we have suppressed the (r, 9) dependence of F. 
Using 



(2.11a) 
(2.11b) 



£{p,a) 

£(p, -a) 
£(a,ip) 



1 /7T 

-- sin — 

2 V2 

1 . / 7T 

-- sm 

2 V 2 

- sm — 
\2aJ 



= - 2 R(p\ 
\r{p), 



as well as the boundary conditions (2.7), equation (2.11a) becomes the following equation: 



q(r,e) 



1 

~2^ 



9+(p) 



In 



P 



tt6»\ 1 



1 



+ 



2vr 



In 



i ii 2 (r) - cos 2 ig) + - A R\p) - X - sin I ^ 1 i?(p)ft(r) 

li? 2 (r) - cos 2 (g) + ^ 2 ( P ) + 1 sin (g) i*(p)i?(r 

/1T(p\ . ( ir8 



dp 



dp 



2 (^) +si " 2 (S)- si HS si Hi |fl(r) 



dyj. (2.12) 



By employing the identity 



1 

1 = 4 



r\ 2^ 
a J 



equation (2.12) becomes equation (2.8a). 



□ 
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Remark 2.1. (Asymptotics) The large r asymptotics of the solution q(r, 9) defined by (2.8a) is given 
by 



q(r,9) 



1 

2^ 



In 



a \a 



S + l »>n(* )( r ) -S + O 

7T \2a J \a. 



where 



and 



5 



s = J [g-(p) - 5+0)] y + a J 9(v)dv 



[g + ( P ) + g-(p)]R(p)j ~ 2« f oMshi (g) dtp. 



oo, (2.13) 

(2.14) 
(2.15) 



Indeed, each of the brackets of the RHS of ( |2.8a[ ) can be written in the form 

'R 2 (r) 



F(R(r),T) =ln 



TR(r) + O(l) 



where T is given, respectively, for each bracket by 



(II i^- 



1 /tt6 
- sin — 

2 I 2a 



R(p), sin 



oo, 



\2aJ \2a) 



As r — > oo, we find 



F(R(r),T) =ln 



R 2 (r) 



1-4T[-) 2a +0 
a 



In ( +21n[i?(r)]-4T(^) J " +0 



InlM + 'ln^ 
4 / a Va 



; / l' r ) 2a +0 

.a 



r — )• oo. 



Using the estimate (2.16), equation (2.8a) implies (2.13). 



(2.16) 



Remark 2.2. (A particular example) Suppose that we choose the functions {g+(p),g-(p),g(ip)} as 
follows: 

a 

Using these functions in the definition ( |2.14[ ) of 5, we find 



g±(p) = {pe ±ia )HP e±ia )> 9(v) = —(ae i *) s u(ae i v), s€ 



(2.17) 



S(s) 



z s 1 u(z)dz, z G C. 



(2.18) 



Indeed, recall that H involves three integrals; making in the first, second and third integrals of the 



RHS of (2.18) the substitutions, respectively, 



z = pe 



z = pe 



z = ae 



up 



we find that the RHS of (2.18) equals — S. 
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Remark 2.3. (The Riemann function) Suppose we chose the functions {g+(p), g~(p), g(<p)} by (2.17), 
where 

u{z) = (2.19) 



e~ z - 1 

Then the function S(s) is proportional to the Riemann zeta function, namely 



S(s) 



2mC,{s) 



z E 



(2.20) 



r(l-a)' 

Hence, the Riemann hypothesis is valid iff there does not exist a solution of the Neumann boundary 
value problem defined in Theorem 2.2 with the functions {g + (p), g~(p), g(y>)} defined by equations 
(2.17) and (2.19), which is bounded as r — > oo. 

The Dirichlet Boundary Values 



By evaluating the RHS of equation (2.8a) at r = a and at 9 = ±a we find the following expressions 
for the Dirichlet boundary values: 



q(a, 



+ 



9+(p) 



In 



g(if) In 



sin 



\2a) 



\2a J 

sm — 
\2a) 



dp + 



7T 



9-{p) 



In 



R(p) + sin 



V2a 



dtp, 



-a < a < a, 



(r, ±a) 



+ 



9+(p) In 
g(tp) In 



R(r) R(p) 
2 2 



do 1 

— + - 

p IT 



9~{p) In 



R(r) R{p) 
2 2 



dp 
P 



R(r) 



=F sm 



/mp 
\2a 



dp, 



(2.21) 



(2.22) 



dcp, 



a < r < oo. 



Adding and substracing equations (2.22) + and (2.22) , we find the following equations: 

R 2 (p) R 2 {r) 



l(r, -a) + q(r,a) 



+ 



7T 



9-{p) -ff+(p)]ln 
R 2 (r) 



dp 
P 



9(f) in 



sm ■ 



2a) 



dip 



(2.23) 



and 



q(r,-a)-q(r,a)=- I [g~(p) + g+(p)) In 



7T 



+ 



R(r)-R(p) 



7T 



g(p) In 



R(r) + R(p) 
i?(r) + 2sin(^) 



dp 
P 



R{r) -2 sin (^) 



(2.24) 



In what follows we assume that the given functions g±(r) and <?(#) satisfy the constraint S = 0. 
Let c be a constant, then 



In 



R(r) 



+ c 



G(r) + 0{r- T / a ), r^oo 



where G(r) is defined by 



G(r) 



a \a 



In 4, r > 0. 



(2.25) 
(2.26) 
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Each of the two logarithmic terms in the RHS of equation (2.23) grows logarithmically as r — >• oo, 



however the condition S = 0, implies that these two terms cancel. Indeed, using the fact that 



G(r)S = 0, equation (2.23) can be rewritten in the form 



1 f°° 

q(r,-a) + q{r,a) =- [g-(p) - g + (p)\ 

* J a 

R 2 {r) 



In 



R 2 (r) R 2 {p) 



G(r) 



dp 
P 



a 

+ - 

7T 



In 



. 7TC9 \- 

sin — 
2a J 



G(r) 



dip. 



(2.27) 



Equations (2.24) and (2.27) are the basic equations needed for the derivation of certain integral 
identities. 

3 The Global Relation and Certain Integral Identities 

The global relations for the Laplace equation in the domain D are the following two equations: 



/ ■ 

JdD 



e ±ik8 r k 



dv 

(—q$ ± ikq) h (rq r — kq)d9 

r 



0, Rek< 



TT 

2a' 



(3.1) 



The above restriction in k is the consequence of the large r behaviour of q, namely of the estimate 

q = 0(r~2^\, T — y oo. 



If the domain D is defined by (2.6), then equations (3.1 1 1 * 1 become 

~e^ iak q(r, -a) - e ±iak q(r, a) 



ret roo 

k{a k e ±ike q(a,9)d6±i / r k 



dr 



a 



k+l 



e ±ike g{d)dd+ I r k e™'j_ (r) - e ±tak g + (r) 



dr 
r 



(3.2) 



Adding and subtracting (3.2^ we obtain the following basic equations: 



k\a k cos(k9)q(a, B)d9 + sin(afc) / r k [q(r, -a) + q(r, a)] 



dr 



,k+l 



cos(k6)g(9)d9 + cos(afc) / r k [g_(r) - g+(r)] 



, dr 



(3.3) 



and 



k\a k sin(k6)q(a,6)d6 + cos(ak) / r k [q(r, -a) - q(r, a)} 



dr 



k+l / sm(k9)g(6)d6 - sm(ak) / r k \g_(r) + g + (r)] 



dr 



(3.4) 



(2.27), we find the following equations, which are valid for Re A; < ^: 



Replacing in the RHS of the relations (3.3) and ( |3.4[ ) q(r,—a) ± q(r, a) by the RHS of (2.24) and 

(3.5) 



a / ai^Fiip^^d^p 



dp 



[g-(p)-g + (p)]F 2 (p,k)^ = 

P 



7 



and 



'I I 0(r)IMr-^'^r + i [9- (fi) + 9+(p)} F*(p, *) ~ = 0, 



(3.6) 



where the functions {Fj}f are defined as follows: 

/a 
cos(k6) In 
-a 

+ k sin(a£;) 
— 7r cos(k(p) 



/7T(p\ I TT0\ 

sin — + sin — 
\2aJ \2aJ 



<:W 



a 



In 



R 2 (r) 



2a J 



G(r) 



-a < tp < a, Re k < 



dr 
r 

TT 

2a' 



(3.7) 



/a 
cos(A;0) In 
-Q 

+ k sm(ak) 



2 V 2 "/ 







f©' 


In 







d0 

i? 2 (r) i? 2 (p) 



dr 
r 



tt cos(ak) ( - 



a < p < oo, Re k < 



TT 

2a' 



(3.8) 



F 3 ((p,k) = -k sin(&0)ln 

J —a 

+ k cos(ak) 
— tt sin(/c(/?) 



. /VTG9\ . f TT0\ 

sm — + sm — 
V2a/ \2a/ 



d6 



-) In 

a 



i?(r) +2 sin ^ 



2a, 



#(r) 



dr 
r 



-a < tp < a, Rek < 



TT 



F A (p,k) =k / sin(/c(9)ln 



2 V 2a / 



7r sin(a/c) 



d6 + k cos(ak) 



a < p < oo, Re A; < 



In 



7T 

2a"' 



2a' 
R(r) - R{p) 



(3.9) 



i?(r) + R(p) 



dr 
r 



(3.10) 



and principal value integrals are assumed if needed. 

The function g~(p) — g+{p) is related to g(ip) through the equation 5 = 0, whereas the function 



g~(p) + g+(p) is independent of g{<p). Hence equations (3.5) and (3.6) imply the following basic 
identities: 

F 1 (<p,k) = f(k), F 2 (p,k) = f(k) (3.11) 

and 

F 3 (tp,k) = 0, F 4 (p,fc) = 0, (3.12) 

where 



a < p < oo, —a < tp < a, Rek < 



TT 

2~a~ 



and f(k) is some function of k. 

The above equations will be verified explicitly in the next section, where it will also be shown 
that the function f(k) is given by 



f(k) = — sin(a/c) 



2 In 2 + 



TT 

ak 



(3.13) 
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Remark 3.1. (Re A; < 0) Using the estimates 



and 



In 



In 



R 2 (r) 



sin — 
2a J 



R 2 (r) R 2 (p) 



a 



G{r) = Q(r 



G(r) = 0{r- 7T/a ), 0<p<oo, r -> oo, 



7T 



it follows that the first two integrals in the RHS of two equation in (3.7) involving the above terms 
are well defined for Re A; < 7r/2a. If k satisfies the stronger restriction Re A; < 0, then it is not 
necessary to subtract the term involving G(r); actually in this case we find 



x / r \ k 



G(r 



dr 



f(k), Rek<0. 



(3.14) 



7T 



k sin(afc) 

Remark 3.2. (Re k < 0) If Re k < 0, then 

Fifak) =F 2 {p,k) = F 3 (ip,k) = F A (p,k) = 0, \ip\ < '- , 0<p<oo, ReA;<0, (3.15) 

where Fi and F2 denote the expression obtained form F\ and F2 by neglecting the terms involving 
G(r). Equations (3.15) can be obtained directly as follows: the function 

(3.16) 



(3.17) 



q(r, 9) = Re (z k j = r k cos(k8) , Re k < 0, 
is a solution of the Neumann boundary value problem with 

9±{ r ) = Tkr k sin(fca), g(6) = /ca fc_1 cos(k9). 
Substituting q(r, a) = r fc cos(fca) and the expressions ( |3.17 ) in equation ( |2.22 ) + , we find 



r k cos(ka) 



k s'm(ka) 



7T 



p k In 



R 2 {p) R 2 {r) 



+ 



ka 



k rot 



7T 



cos(kip) In 



4 

R(r) 



sin 



dp 
P 

/TUp\ 

\2a) 



dip, Re k < 0. 



(3.18) 



Replacing ip with — y and dividing by a k /tt, we find F2 = 0. 

Similarly, substituting q(a,9) = a k cos(k9) and the expressions (3.17) in equation (2.21), we find 
F 1 = 0. 

The function 

q(r, 9) = Im (V) = r k sin(k0), Rek<0, 
is also a solution of the Neumann boundary value problem with 

g±(r) = ±kr k sin(ka), g(9) = ka k ~ l cos(k9). 

The above solution implies F3 = F^ = 0. 
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Remark 3.3. The first integral in the definition of Fi involves a singularity at 9 = —<p. Making the 
change of variables 9 = — ip + x, the relevant singularity is mapped at x = and the associated 
integral is ln|x|. This singularity can be handled by Cauchy principal value integrals. In this case 
the contribution of this singularity vanishes. Indeed, this contribution involves 



lim 



+ / I In \x\dx = 2 lim 



In xdx 



and using integration by parts it follows that the contribution from x = e vanishes. 

If instead of the principal value integral we use the limit from above and below, the relevant 
contribution is 



limi I \u{e)e id d9 



which clearly does not exist. 



4 Verification of the Four Identities 



In order to simplify the functions {Fj}f defined by (3.7)-(3.10), we introduce the change of variables 



r _2a 

- = X T , 

a 



P _2q 

- = y 

a 

u = e 2t » , v = e 2 « . 



Then {Fj}f can be written as 



—J [h 3 {n, v) + h 3 ( K , -v)} - -(v K + v~ K ), 

fTTK\ / UK 

F 2 = -Kh 6 (K,y) +KSm ^ — j [h A {K,y) + h 5 {K,y)\ - Try cos ^ — 



—J [h 3 (K, v) - h 3 ( K , -v)] + 2~^ vK ~ V ~ K "> 

F A 

where the functions {hj} 1 are defined by 



-Kh 7 (n, y)-K cos (J^-j [h^K, y) - h 5 (K, y)) - iry K sin (™ 



hi(K,v) 
h 2 (K,v) 
h 3 (K,v) 
hi(K,y) 
hs(n,y) 
he{K,y) 
h r (n,y) 



71 



U K 1 + U K 1 

2i 

U K-1 _ U -K-1 



v — v 1 u — u 1 



2i 



+ 



2i 



71 
1 



v — V 1 u — u 1 

+ 



2/ 



du, 
du, 



„K-1 



x + x 1 v — v 1 



2 

X + x~ 



2i 

y + y~ 



2 

x + x~ 



+ 



2 

y + y~ 



U K 1 + u K 1 



2 2 



2/ 

U K-1 _ U -K-1 



71 



2 



+ 



2i 

+ In x + In 2 
+ In x + In 2 
+ In x + In 2 
du. 



dx, 
dx, 
dx, 



u — u 



+ 



2i 
u — u 



2i 



du, 



(4.1) 

(4.2 
(4.3 
(4.4 
(4.5 

(4.6 
(4.7 
(4.8 
(4.9 
(4.10 
(4.11 
(4.12 
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and 71 is depicted in figure |4.1| 

In order to compute {hj}\, the following lemma will be useful. 

Lemma 4.1. Let 7 be a smooth curve from z\ to z 2 in the z— complex plane which passes through 
zq . Then 



I = PV 



Z - Z 



-dz 



±—[z p+1 2 F 1 (l,p+l; P + 2 ] 
■ i-iZn 



(p+l)*b 



Zi, 



- z p ^ 2 Fx{l,p+ l;p + 2; ^)] - niz p , 



Rep > -1. 



(4.13) 



// i/ie singularity lies off the path 7, then the last term, in (4-13) is absent. 
Proof. Letting x = z/zq, we find 



I = z p Q PV 



70 



x - 1 



-dx, 



where 70 runs from zi/zq to z 2 /zq, avoids x = and passes through x = 1. We recall that 



p + 

Letting px = v, we find 
so that 



2 Fi(l,p+l;p + 2;a;) = / -^—dp, Rep > -1. 
-f 1 J 1- px 

1 /-x „p 

j 2 Fi(l,p+l;p + 2;x) = y - 



p + 

p 



1 d 



x — 1 p + 1 dx 



[x p+1 2 Fi(l,p + l;p + 2;x)], 



(4.14) 



which immediately gives (4.13), the principal part at zq being included by subtracting half the residue 
at this point. □ 



The hypergeometric functions in (4.13) are related to the Lerch transcendent, a generalization of 
the Riemann Zeta function 

2 F 1 (l,p;p+ l;z) =p$(z,l,p); (4.15) 
the Lerch transcendent is analytic in the z-plane for fixed p 7^ 1. 
Proposition 4.2. For Re k > —1, 



F\ (f,k) = F 2 (p, k) = - sin(afc) 



2 In 2 + 



7T 

ak 



(4.16) 



Proof. If k = (i.e. k = 0), then Fx = F 2 = — tt, which agrees with (4.16) in the limit k — > 0. Thus, 
we consider the case of k ^ 0. Using integration by parts and a partial fraction decomposition, we 
find 



hx(K,v) =- sin 

K \ 2 



hi 



+ 



1 



2in(v + v 1 ) 



v — v 
PV 



-1 



+ 1 



+ In 



v — v 
~27 



-1 



1 



Ut — u 



-)(u + 



71 



1 



1 



u + u u — f 



1 dit, 



(4.17) 



where the contour 71 is depicted in figure |4.1| and we have used 



(i) K -H) K = 2isin 



11 



Im u 




Figure 4.1: The oriented contour 71. 



Using the change of variables u — > 1/u for the term involving u K , the above principal value integral 
becomes 



2in{v 



^PV [ u^u + u- 1 ) (— 1 ^— + ^—) 

+ v 1 j J yi \u + v u — v 1 u — v U + V~ L J 



du. 



(4.18) 



Each integral in (4.18) can be computed by using lemma 4.1 Noting that —a < <p < a, the definition 
of v (see (4.1 )) implies that v and v _1 lie on the contour 71. The relevant residue contributions yield 
the term 



7T 



(v K + v~ K ). 



(4.19) 



Hence, we find 

-PV [ u K (u + u- 1 ) (— 



2in{v 



+ v u — v 



-1 



1 1 

+ 



u — v U + V 



-1 



du 



1 /7Tfv 

— tt sin — 

k(v + v- 1 ) V 2 



1 _ 1 _ 

-F(k + 1; -i/u) + -F(k + 1; i/v) + i>F(re + 1; -iv) 

V V 



+vF(k + 1; iv) F(k — 1; —i/v) F(re — 1; i/u) — vF(k — 1; — iv) — vF(k — 1; 



7T 

2^ 



(4.20) 



where 



F(k;u ) 



1 



K+ 1 



2^1(1, re + 1; re + 2; u ) 



(4.21) 



Regarding hs(K,v), using integration by parts and a partial fraction decomposition, we find 



1 


jln 


f — u 1 






K 




2i 






/■1 







1 +ln2 

re 



X + IV x — IV 



dx > . 



(4.22) 
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i 

iy 



-iy 



-i • 

Figure 4.2: The deformed contour 72. 



Each integral of the second line in (4.22) can be computed using lemma 4.1 The definition of v 



implies that the integrand in (4.22) does not have any singularities on [0, 1]. Hence, we find 

1 



h 3 (K,v) 



1 



jln 


v — v 1 






2i 



+ ln2 



1 



v + v 



-1 



-F(k + 1; i/v) - -F(k - 1; i/v) 

V V 



+vF(k + 1; — iv) — vF(k — 1; —iv) 



(4.23) 



Substituting (4.20) into (4.17) and using (4.23), we find that the terms involving the hypergeo- 
metric functions and the logarithmic functions cancel. Hence, recalling k = — — k, we obtain (4.16) 
for Fi. 

Regarding F2, we first deform the contour 71 to 72, which is the segment (— i, i) of the imaginary 

; integration b 

(y + y _1 ) + l 



axis, see figure |4.2[ Then, using integration by parts and a partial fraction decomposition, we find 

1 



1 . /7TK\ 



h(K,y) =-sin j 



In 



+ In 



(y + y- 1 )-! 



+ ^7^/ 



1 1 

— + r 



iy 11 + iy 1 u — iy ' u — iy 1 



du. 
(4.24) 



Using the fact that the integrand of the above principal value integral has poles at x = ±iy and 



employing lemma 4.1, we find that the principal value integral of the RHS of (4.24) is given by 
1 



n(y - y 
1 



TTK 

. . sin ! 

V 2 



( 



-F(k + 1, -1/y) - yFlyK + 1, -y) + -F(k + 1, 1/y) - yF( K + 1, y) 

y y 



1 



F(k - 1, -1/y) + yF( K - 1, -y) - -F{k - 1, 1/y) + - 1, y) 

y y 



IT fKK\ 



(4.25) 
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For h±(K,y) and h^{K,y), using the fact that the integrands have a pole at x = y, we find 



h^K, v) 



In 



y + y 



-i 



1 , 1 

- + ln2 + 



-yF(K + l\y) + yF{K-l;y) 



y-y 

+ iwy 



-F( K + l;l/y)--F(K-l;l/y) 

y y 



(4.26) 



and 



/i 5 (k,i>) = - <! In 

K 



y + y 



-i 



+ i 



l l 

- +ln2 + 



y-y 



-i 



-F( K + l;-l/y)--F(K-l;-l/y) 

y y 



-yF{K + l;-y)+yF(K-l;-y) 



(4.27) 



Combining (|4.24|), (|4.26|) and (|4.27|) with (|4.25|), we find that the terms involving the hypergeo- 

□ 



metric functions and the logarithmic functions cancel and hence we find (4.16) for Fi 
Proposition 4.3. For Re« > — 1, 

F 3 ((p,k) = F 4 (p,k)=0, 



(4.28) 



Proof. If k = 0, it is obvious that F3 = F4 = 0. Thus, we consider the case Using integration 

by parts and a partial fraction decomposition, ^(k, v) can be written as 



h 2 (K,v) =- cos (J^j 



hi 



v — V 



2i 



+ 1 



In 



v — v 



-1 



2i 



1 



+ 7T^ tv PV [ u^iu + u- 1 ) 



1 



1 



1 



u + v u — v 1 u — V u + v 



T du, (4.29) 



where we have used 



(i) K + (-i) K = 2cos 



TTK 



The principal value integral in ( 4.29| ) can be evaluated by lemma 4.1 
1 



/TTK 

—, r- cos — 

^(v + v- 1 ) \ 2 



1 _ 1 . 

-F(k + 1; -i/v) F(k + 1; i/v) - vF(n + 1; -iv) + vF(k + 1; iv) 

V V 



1 _ 1 „ - _ 

--F(k - 1; -i/v) + -F(k - 1; i/v) + vF(k - 1; -iv) - vF(k - 1; iv) 

V V 



IT 

2in 



(v K -v- K ). 



(4.30) 



Substituting (4.30) into (4.29) and combining the resulting expression with (4.23), we find F$ = 



Regarding F4, deforming the contour 71 to 72 and then using integration by parts we find 



1 




In 


— COS 


:f) 


K 







(y + y x ) + i 



In 



-(y + y" 1 )-! 



nry 



+ -A vtPV [ u K (u + u- l )(— 



1 1 

+ 



iy u + iy 1 u — iy u — iy 1 



du. 
(4.31) 
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Using the fact that the integrand of the above principal value integral has poles at x = ±iy and 
employing lemma 4.1, we find that the principal value integral of the RHS of (4.31) is given by 

1 



/ TTK 

— ; TV COS 

Kiy-y- 1 ) \2 



1 F(k + 1, -1/y) - yF( K + 1, -y) - 1 F(k + 1, l/y) + yF( K + 1, y) 

y y 



lp( K - 1, -l/y) + y P( K - 1, -y) + -F(k - 1, l/y) - yF(K - 1, y) 

y y 



IT / TTK 

— y sin 

K 



Using (4.26) and (4.27) together with (4.32), we find F 4 = 0. 



(4.32) 

□ 



5 Functional Identities 

The functions /13, /15, and /17 do not involve principal value integrals, which always yields a term 
involving i, thus for k real, it should be possible to express these functions in terms of real functions. 
We have succeeded in doing this for the last three functions but not for /13. The relevant expressions 
are 

1 



h 5 (k,y) 



he(k,y) 



h 7 (k,y) 



Uln ^(y + y-^ + l 



k 
+ 



k 



+ ln2 



1 



Hy- 2 - 1) 
1 



[y~ 2 2 F 1 (l, k; k + 1; -y- 1 ) - 2 F l (l, k;k + l; -y)) 



(fc + 2)(y-2-l) 

(y-y- 1 ) 2 In 



[y- 2 2 i ? i(l,fc + 2;A; + 3;-y- 1 )- 2 *i(l, k + 2; k + 3; -y)] \, (5.1) 



sin(f ) 
k(y - y 



-l^2 



.(y-y^ 1 ) 



fc-2 



3-^2 



1 1 l . 2 - k k - 4 

2' ' ' 2'2' (y-y- 1 ) 2 



+ 



1 



/c + 2 
2 



3-^2 



1 k k 4 

-,1,1;--, 2+-;--— — ^ 



(5.2) 



fc 



cos 



+ 



y + y~ 



irk\ f 1 
Tj\2 
1 



In 



2/ + y 



-1 



(y-y x ) 2 
1 



k + l 



3F2 



2' ' ' 2 



In 



1 - k 3 + k 



y + y 



-1 



+ 1 



2 ' {y-y^Y 



13-fcl + Jfe 4 
//-l 3 2 V 2' ' ; "2"' ~ (y - y- 1 ) 2 



Refc>-1, lny<0. (5.3) 



In addition, for k real, the functions h\ and h 2 possess the following more complicated real forms: 



I) ] ( /,-. /' ) =— -(2 sin (^^j ^ n 

(v + v' 1 ) 2 



and 



h 2 (k,v) 



k 



v + v 



+ 



IT 



3/2 



r<2i 

U 33 



firk\ 
cos \Y ) 



In 



1 2 - 
1 1 

f — v~ 



v + v *) 2 

Jfc k 
2 J 



U 33 



(v + t^ 1 ) 2 



1 -I 2 + | 



1 1 



(5.4) 



2i 



In 



v — v 



-1 



+ 1 



2^ 2 ( 



v — v 



+ G33 



i(v + v 1 ) 
(v + -u -1 ) 2 



G 



1\2 



(v + W _i ) 



21 

l-fc 3+fc 
2 2 

1 1 2 



1 



3-fc 1+fc 



2 

1 1 



(5.5) 
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where G™ n denotes the Meijer G-function. 



The basic identities (3.11) and (3.12) yield a plethora of novel identities involving the hypergeo- 



metric and related functions. 



Example 5.1. Replacing in the definition (4.5) of F4, h-j by (5.3), h± by (4.26) and h*> by (5.1), the 
identity F4 = yields: 



+ 2,-1)2 

,,k 

i + tan 



.«|i,i,ii 1 -* 8 + * 



Try" 



y + y 1 

1 



2 
irk 

T 



2 ' 2 ' (y-y-^' +{k ^ k) 



+ 



V 



2 F 1 (1, k; 1 + k;y) - 2 F 1 (1, fc; 1 + k; -y) 



k + 2 



y 2 — y 2 \ k 
2 1- 1 ( .1 ■ A" I- 2: A- + 3; y) - 2 F 1 (1, k + 2: A- + 3: - // ) 

Refc>-1, lny<0. (5.6) 



The i appears above as a consequence of the fact that the hypergeometric function acquires an 
imaginary part when continued to real arguments exceeding unity. 



Letting k = 2 with z = 2/(y — y 1 ) equation (5.6) yields the following novel identity: 



^2 ^,1,1; 



1 5 
2'2 ; 



which, by analytic continuation, is valid for all z. 



3 f 4 + 3z 2 sinlT 1 z 



z 2 Wl+^ 



(5.7) 



For k = 1/2, equation (5.6) yields the curious identity 
4 



{y-y x ) 2 



3F2I 3 



\ 5 ' ' A ' >* 5 



3(1 + f)^ 3 / 2 ^ 
+ 6- 



y 2 + l 



y 2 + l 



f 1 3 5 4 
1" 1 ( <y - .^Z- 1 ) 2 ; + 3 2 U' 5 4' 4 5 " (y - 

//I tiinJi — - — tan -1 ] + ( tanh -1 ,/y — tan -1 ,/y ] 

V y/y VVJ V / 



/ 1 tanh 1 — — 



(5.8) 



for < y < 1, where since the left hand side depends only on 4/ (y — y -1 ) 2 , so must the right hand 
side, which is not at all obvious. 

Similarly, for k = 4, we find the novel identity 



/ I 1 1 I-- 3 7 - 4 

1 2' 2' 2' (y-y 1 ) 2 

5(y-y- 1 ) 2 (y 4 + y 2 + l 



y 

y 2 < 1. 



2(2/2 + 1) 



sinh 1 



2y 



y 2 + l 1 fl + v 

+ o , o TV m 



l-y 2 J 2y 2 {y 2 + l) \l - y 



1 1 



(5.9) 



Example 5.2. Replacing in the definition (4.3) of F2, h% by (5.2), /14 by (4.26) and h§ by (5.1), the 
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identity F2 = f(k) yields the following identity: 



{k + 2 )(y-y-iy 



! ( J, 1,1;~ 2 + ifc; 



27ry fc z - cot 



irk 



+ 



2' 2 
2y- x f 1 



+ (fc ->■ -fc) 



-1 



(y-y- 1 ) 2 , 

2 Fi(l, fc; fc+ 1; 2/" 1 ) + 2^1(1, k;k + l; -y~ l ) 



k + 2 



2 Fi(l, fc + 2; fc + 3; y" 1 ) + 2 Fi(l, 2 + fc; 3 + fc; -y' 1 ] 

Rek > -1, 



lny| < 00. (5.10) 



Replacing in the definition ( 4.4 ) of F4, /i 2 by ( 5.5 ) , we find the novel identity, valid for Re k > — 1 
and all real v, 



G 



(v + v~ r ) 2 



1 



l+k 3-k 



J I I )+(k^-k) 



(v + v- 1 ) 



m 1 l' i (v — v 1 ) I vr 
1 



^ 2 <2 firk 
— cos 



1 2v 

— — Re[— f 2 F 1 (1, 2 + k;3 + k;iv)-(v^ v~ l ) 

k + 2 v + v 1 



fc + 1 



Re 



v 2 - 1 



£(« — f 1 ) 



2 F X (1, fc + 1; fc + 2; iv' 1 ) - (v -> u" 1 ) 



In particular, for fc = and all real 6, 



G§£ ( cos 2 6 



1 1 2 
x 2 2 

1 1 i 



2 , , , / cos 6 

-—7^ cos cot in 

vr 3 / 2 Vl + sin6 



(5.11) 



(5.12) 



The above identities appear new. Thus, it seems that by employing the global relation to the 
solution of certain boundary value problems, it is possible to construct new formulas in the area 
of special functions, although it is not clear how the form of these identities can be predicted in 
advance. 
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